Spin dynamics of molecular nanomagnets fully unraveled by four-dimensional inelastic 

neutron scattering 
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Molecular nanomagnets are among the first examples of spin systems of finite size and have 
been test-beds for addressing a range of elusive but important phenomena in quantum dynamics. 
In fact, for short-enough timescales the spin wavefunctions evolve coherently according to the an 
appropriate cluster spin-Hamiltonian, whose structure can be tailored at the synthetic level to meet 
specific requirements. Unfortunately, to this point it has been impossible to determine the spin 
dynamics directly. If the molecule is sufficiently simple, the spin motion can be indirectly assessed 
by an approximate model Hamiltonian fitted to experimental measurements of various types. 
Here we show that recently-developed instrumentation yields the four-dimensional inelastic-neutron 
scattering function 5*(Q, E) in vast portions of reciprocal space and enables the spin dynamics to 
be determined with no need of any model Hamiltonian. We exploit the Cr8 antiferromagnetic ring 
as a benchmark to demonstrate the potential of this new approach. For the first time we extract 
a model-free picture of the quantum dynamics of a molecular nanomagnet. This allows us, for 
example, to examine how a quantum fluctuation propagates along the ring and to directly test the 
degree of validity of the Neel- vector-tunneling description of the spin dynamics. 
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Mesoscopic systems can exhibit typical quantum dy- 
namical phenomena, for instance by being able to tunnel 
through an energy barrier or by displaying long-lived 
coherent oscillations associated with superposition of 
states. This has attracted considerable interest for 
addressing fundamental issues and for the possible appli- 
cations in quantum-information processing. Molecular 
nanomagnets (MNMs) are spin clusters where the topol- 
ogy of magnetic interactions can be tailored precisely at 
the synthetic level. They are metal-organic molecules 
containing a small number of magnetic ions whose spins 
are strongly coupled by exchange interactions. Shells 
of organic ligands provide magnetic separation between 
adjacent magnetic cores, which behave as identical and 
independent zero-dimensional units [l[. The magnetic 
dynamics are characterized by strong quantum fluc- 
tuations and this makes MNMs of great interest in 
quantum magnetism as model systems to investigate 
a range of phenomena, such as quantum-tunnelling 
of the magnetization 043 j Neel- vector tunnelling 
(NVT)0, Q, quantum information processing_0J3 , 
quantum entanglement [lol - [l3j or decoherence [l4 - |l6f . 
Besides their fundamental interest, MNMs are also the 
focus of intense research for the potential technological 
applications as classical or quantum bitsp], 0"S [13] an d 
as magnetocaloric refrigerants flij. 

A crucial aspect of the research on MNMs is the 
understanding of their low-temperature spin dynamics, 
especially of those aspects which are a direct manifesta- 
tion of quantum mechanics like the tunneling of the Neel 
vector in antiferromagnetic rings. The most powerful 
technique to investigate the spin dynamics is inelastic 



neutron scattering (INS). INS measurements on MNMs 
are normally performed on polycrystallinc samples. For 
these the cross-section is averaged over all molecular 
orientations, and only depends on the modulus of 
momentum transfer (Q). The resulting INS spectrum is 
composed of peaks whose energies yield the eigenvalues 
and whose areas contain information on the correspond- 
ing wavefunctions. The dynamics is then extrapolated 
by fitting experimental spectra to a suitable spin Hamil- 
tonian H and by performing calculations within the 
framework set by this model. However, this indirect 
approach does not work in many interesting cases when 
the Hilbert space dimension is too large to diagonalise H 
or when a large number of model parameters makes the 
fit non-univocal. Even by using single-crystal samples, 
only limited and partially-integrated information on 
the Q-dcpendcncc of the scattering function S(Q,u>) 
can be obtained by traditional time-of-flight neutron 
spectrometers, i.e., carrying unitary detectors fitted 
on Debye-Scherrer rings. The implementation of large 
arrays of position sensitive detectors in cold-neutron 
time of flight spectrometers fill l2fj| . together with the 
advances in software [21], has recently opened un- 



precedented possibilities in inelastic neutron scattering 
experiments on molecular nanomagnets allowing to 
determine the four-dimensional scattering function 
S(Q,lu) in a vast portion of the reciprocal space. This 
provides a much more selective characterization of the 
MNMs when different candidate models can only be 
discriminated by the vectorial Q-depcndence of S(Q,ui) 

Here we show that the measurement of S(Q,lu) over a 
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FIG. 2: Constant energy plots of the neutron scatter- 
ing intensity for the p=l excitation. Cuts are shown in 
the Qx'-Qy', Qz'-Qy' and Qz>-Qx> planes, where the primed 
reference frame is defined in Fig. SI (Supplementary informa- 
tion). The measurement was carried out with a 5 A incident 
neutron wavelength, and a sample temperature of 1.5 K. 



I. THE Cr 8 PROTOTYPE 
ANTIFERROMAGNETIC RING 



FIG. 1: Magnetic energy spectrum of Crg and zero- 
temperature inelastic neutron scattering transitions. 

(a): Low- lying energy multiplets as a function of their to- 
tal spin for the isotropic exchange Hamiltonian of Cts (Eq. 
|T|) with D = 0). Arrows indicate the three transitions 
seen by INS at zero temperature. All other transitions have 
negligible cross-section. Blue and red symbols indicate L- 
and S-band states respectively. The inset shows the core of 
Crg (C8oCrsFgDi44 032 ; green=Cr, yellow=F, red=0, dark 
grey=C, D omitted), (b): measured low-T INS spectra for 
Crs with an incident neutron wavelength A = 3.1 A. The la- 
bels indicate the three peaks corresponding to the transitions 
reported in (a). The p = 1 transition is partially hidden by 
the elastic signal. The inset reports higher-resolution mea- 
surements with A = 5A showing the p — 1 transition split by 
magnetic anisotropy. 



large solid angle of wavevectors in high-quality single- 
crystals of MNMs permits us to directly determine the 
full pattern of real-space dynamical two-spin correlations 
without assuming any underlying model Hamiltonian. 
These correlations are the key quantities characterizing 
the magnetic dynamics. For instance, they determine the 
linear response of the system to an arbitrary magnetic 
field varying in space and time. We exploit Crs, a proto- 
type antiferromagnetic (AF) ring (Fig. [T]and Methods), 
as a benchmark to demonstrate the potential of this new 
approach for the determination of the spin-dynamics 
in MNMs. The so-obtained correlation functions allow 
a model-free picture of the quantum dynamics of the 
ring. For example, the way a quantum fluctuation 
propagates along the ring or the degree of validity of 
a Neel vector tunneling description [El, [23[ are determined. 



The Cr§ molecule contains eight Cr 3+ ions (s = 3/2) 
forming a nearly regular octagon [24|, [2f| . It is one of the 
most studied MNMs both for its intrinsic interest and 
for being the precursor of a large family of heterometal- 
lic molecules and supramolecular complexes which are 
intensively studied for a number of fundamental and ap- 
plicative issues 0, [H, [ID, l26h 28| . The spin-Hamiltonian 
of Crs contains a dominating nearest-neighbor Hciscn- 
bcrg AF exchange and small axial anisotropic terms: 

8 8 

H = J^s( l ).s(z + l) + J D^ S 2 (l) (1) 

i=l i=l 

where s(9) = s(l), J = 1.46 mcV, D = -0.038 meV 
and the anisotropy axis z is perpendicular to the 
ring plane[H HI- The ground state of Eq. CO is 
a nonmagnetic singlet (total spin 5 = 0), and the 
low-lying excited levels are arranged into rotational 
bands (FigfTJ) . The lowest one (called L-band) contains 
the ground state and approximately follows the Lande 
rule E(S) = 2JS(S + l)/iV[J H|. The second set 
of levels belongs to the so-called i?-bands, which are 
also parabolic with respect to S but shifted to higher 
energies. Due to their internal structure, L-band states 
can be excited by neutrons practically only to L- or 
.E-band statesdj. Hence at low temperature, where 
only the S = ground state is populated, three peaks are 
expected and observed in INS spectra[25( (transitions 
marked by arrows in Fig. [1]). Anisotropy produces 
small splittings of the otherwise degenerate 5-multiplcts 
and a tiny second-order mixing of different multiplets. 
For instance, the L-band S = 1 triplet is split into an 
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a 0.66 < E < 0.98 meV 



b 3.8 < E < 4.2 meV 
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FIG. 3: Constant energy plots of the neutron scattering intensity showing all the possible magnetic excitations 

in Crs at 1.5 K. We show the dependency of intensity on two wavevector components Qx-Qy lying in the ring's x-y plane, 
integrated over the full Q z data range, a: Data from transition p = 1, measured with a 5.0A incident neutron wavelength, 
b-c, Data from transitions p = 2 (b) and p — 3 (c), measured with a 3.1 A incident neutron wavelength, d-f, fits to Eq. @ for 
excitations 1 (d), 2 (e) and 3 (f). 



\S = 1,M = 0) singlet and an \S = 1,M = ±1) doublet 
which are resolved by high-resolution INS measurements 
(inset of Fig. [TJd). Splittings of the i?-band triplets are 
smaller and not experimentally resolved. 



II. INELASTIC NEUTRON SCATTERING AND 
DYNAMICAL CORRELATION FUNCTIONS 

The INS spectra in Fig. [Tp are integrated over a wide 
range of momentum transfer Q and hence contain no 
direct information on the spatial structure of wavefunc- 
tions. This is usually indirectly inferred by fitting the 
spectra to a specific model Hamiltonian like Eq. [T] By 
exploiting the new position-sensitive-detectors setup of 
the cold- neutron time of flight spectrometer IN5[l||, we 
have measured the detailed Q-dependence of the scat- 
tering function S(Q,uj) in a vast portion of reciprocal 
space. Fig. [2] and the three upper panels of Fig|3] display 
examples of the Q-dependcnce of S(Q,u>) for the three 
low-temperature INS transitions. Since we would like to 



first address the effects of isotropic exchange, data for the 
p = 1 transition are integrated over the energy range of 
the anisotropic splitting. Anisotropy will be addressed in 
the last part of the paper. Fig|3] shows that the variation 
of S^Q, to) in the Q x -Q y plane is characterized by several 
maxima and minima whose shape and positions depend 
on the specific transitions and reflect the structure of the 
involved wavefunctions. More specifically, for a MNM 
with uniaxial anisotropy and T — > [32] 



a=x,y,z \ ^ ' p d>d'=l 

cos(Q ■ R dd ,)(0\s a {d)\p)(p\ Sa {d')\0)5{cj - E p /h) (2) 

where Fd(Q) is the form factor for the dth ion [§31, |0) 
and \p) are ground and excited eigenstates, respectively, 
E p are eigenenergies and Rdd' are the relative positions 
of the N magnetic ions within a molecule. As is well 
known, this formula can be rewritten in terms of T = 
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FIG. 4: Propagation of a local disturbance deduced from the present INS spectra and in a Neel-vector tunneling 
regime. The frames show the time evolution of (s z (i)) in an eight-spins ring after a delta-pulse perturbation —bs z (l)8(t) applied 
on the red (d = 1) spin. Note that the ground state being a singlet, (s z (i)) = just before t = (frame not shown). Top 
: dynamics deduced from Eq. (5) using the present experimental Fourier coefficients (Table I) and frequencies. The delay 
between two frames is 1.9 10 -13 s, i.e., l/32th of the longest oscillation period 2ir/u)i. Bottom : results for a hypothetical ring 
in a Neel-vector tunneling regime, i.e., D/J = —0.2 corresponding to a tunnel action So ~ 7.6. The delay between two frames 
is l/32th of the tunneling oscillation period 2-kTi/ 'A, A being the tunneling gap. 



dynamical correlation functions 

( Sa (d,t)s a (d\0))=Y,(0\^(d)\ P )(p\s a (d')\0)e-^\ 
p 

(3) 

where lu p = E p /fr. In fact the Fourier coefficients 
c dd'( w p) = (0\s a {d)\p)(p\s a (d')\0) in Eq. © coincide 
with the coefficients in Eq. @. While the values of 
uip are directly read out from the energies of the peaks in 
the INS spectrum, the c dd , (ui p ) can be extracted from the 
data by fitting Eq. ^ to the observed Q-dcpendence of 
each peak. Indeed, for each value of E p these coefficients 
are the only unknown quantities in Eq. (|2|). In order to 
fix the scale-factor in Eq. [5] we have exploited the sum 
rule 

(0| S 2 (d)|0> = S ( S + l) = ^c^K) ( 4 ) 

a,p 

with s(s + 1) = 3.75. 

The Fourier coefficients c dd ,(oj p ) extracted by the fits of 
data integrated over the anisotropy splittings are given in 
Table I. They do not depend significantly on a within ex- 
perimental uncertainty, consistently with the first-order- 
perturbation effect of anisotropy in Crs. The three lower 
panels of Figj3] show calculated Q-dependencies corre- 
sponding to the best- fit Fourier coefficients c" d (u> p ). The 
agreement with experimental maps is excellent. The in- 
dependence of c" d (u> p ), which determines the spatial pat- 
tern of correlations at the various frequencies, strongly 
varies with uj p . This is in line with the behavior expected 
from the model of Eq. Q] (see Table I) . 
The information on the low-T spin dynamics embedded 
in the dynamical correlation functions (Eq. ([3])) can 
be visualized by exploiting the well-known link between 
these correlations and linear response functions [32|. The 
building blocks of such functions are the set of suscepti- 
bilities Xa d which provide the response of s a (d) at 
time t to a delta-pulse perturbation produced by a field 



TABLE I: Fourier components of two-spin dynamical correla- 
tions Ci d (u} p ) — cl d (u} p ) = c\ d {ui p ) extracted from the present 
data (Exp column), and calculated by the isotropic Heisen- 
berg model (Eq. ([T]) with D = 0). The uncertainty on the 
experimental coefficients is of the order of 0.04 for d > 1 and 
of 0.06 for d = 1. Coefficients for d = 6, 7, 8 are equivalent by 
symmetry to d = 4, 3, 2, respectively. 
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b of direction a' applied on spin d' : 

SH(t) = -bs a ,(d')5(t) -+ (s a (d)) = bxtfit) (5) 

with 

xt;f{t) - l j.®{t){K{d,t), Sa ,{d')]) 
2 

= 5 a ,a'x9(<)X) C dd'( w p) sin ( w P*)> ( 6 ) 
p 

& being the step function. The response here is diagonal 
(8a,a') because of axial symmetry. As an example, we 
show in Fig. [4] how a fluctuation on a given site at 
t = 0, (s z (l)) > (produced by a delta-pulse field on 
site 1) propagates in subsequent times along the ring. 
The three-frequencies spectrum of Table I produces a 
wavclike motion of the magnetization, which propagates 
forth and back along the ring with a pattern determined 
by the precise values of Fourier coefficients. After the 
pulse is applied on site 1, waves propagating clockwise 
and counterclockwise interfere constructively at the 
opposite (d = 5) site, thus producing a large value of 
local magnetization (s 2 (5)) after about 6 x 10~ 13 s. 
Such constructive interference is lacking in the case 
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of an odd-numbered AF ring. Here frustration would 
produce a complex time dependence when the two waves 
meet, which could be directly deduced by the present 
technique. 

The measurement of the whole set of low-T Fourier 
coefficients also allows us to extract equal-time cor- 
relation functions (i.e., Eq. ([3]) with t = 0), which 
arc important quantities widely used to characterize 
the ground-state. The advantage here is that one can 
isolate the weak magnetic signal by integrating over the 
inelastic response. Such information would be extremely 
challenging to obtain directly by other techniques 
like diffuse scattering, even with polarization analysis. 
The staggered pattern of equal-time correlations (see 
Supplementary Information, Fig. S2) has an envelope 
slowly decreasing with distance as expected for the one- 
dimensional AF Heisenberg model with s = 3/2, which 
has a quantum critical ground state with power-law 
correlations [34) . 



III. ANISOTROPY 

As discussed above, magnetic anisotropy in Crs is 
small and its main effect is to cause a splitting of the 
S = 1 triplet into a \S = 1,M = 0} singlet and a 
\S = 1,M = ±1) doublet. This splitting is witnessed 
by the presence of two separate peaks at 0.7 and 0.9 
meV in the high-resolution INS spectrum (see the inset 
of Fig. [T}d). We have investigated the effects of mag- 
netic anisotropy on the spin dynamics by measuring the 
Q-dependence of the neutron scattering intensity sepa- 
rately for the two low-energy peaks. The results are 
shown in the two upper panels of Figure S3 (see Sup- 
plementary information). The positions of the maxima 
are the same for both peaks but the relative intensities 
of the various maxima arc different. In particular, the 
two Q x = maxima are the most intense for the 0.7 
meV transition but are the weakest for the 0.9 meV one. 
We have followed the same procedure as before to fit 
separately the Q-dependencies of the two peaks. The 
so-extracted Fourier coefficients are reported in Table II 
and the corresponding calculated maps are displayed in 
the two lower panels of Figure S3. Table II shows that 
c< dd'i. UI p) are non-negligible only for a = z for the first 
peak and for a = x,y for the second peak. This is ex- 
actly the behavior expected for the axial Hamiltonian of 
Eq. {!]), and means that the low- frequency oscillations of 
(s z (d, t)s z (d', 0)} occur at a somewhat smaller frequency 
than those of {s x (d, t)s x (d', 0)) and {s y (d,t)s y (d' ,0)) . We 
find a slightly larger weight of the zz coefficients (first 
peak) with respect to the xx and yy ones (second peak), 
consistent with the asymmetry expected from the model 
of Eq. JT]) in presence of S- mixing [25j. 



TABLE II: Fourier components of two-spin dynamical corre- 
lations c± d (ujp) = c\ d (uj p ) and cl d (u} p ) directly extracted from 
the present data. The uncertainty on the experimental coef- 
ficients is of the order of 0.05. Calculations with the model 
of Eq. Ul yield: alternating ±0.86 for the zz coefficients and 
zero for the xx and yy ones (0.7 meV peak); alternating ±0.73 
for the xx and yy coefficients and zero for the zz ones (0.9 
meV peak). The average of the xx, yy and zz coefficients is 
consistent with the values for the unresolved p = 1 transition 
in Table I. 
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IV. NEEL- VECTOR TUNNELING 

It has been proposed [f| that the low-temperature 
dynamics of antiferromagnetic rings could be charac- 
terized by the coherent tunneling of the Neel vector 
n = Yld(~ l) d s(d) between the +z and —z directions. 
This phenomenon is the AF counterpart of the tunneling 
of the magnetization (QTM) observed in single-molecule 
magnets Interestingly, while for the latter the 

tunneling time is so long that the dynamics becomes 
overdamped due to dissipation into the nuclear-spins 
subsystem, the time scale of NVT in AF rings could 
be much shorter than the decoherence time. Within 
the semiclassical framework, the condition for the 
occurrence of Neel vector tunneling is that the action 
Sq = Nsy/—2D/J be much larger than unity. Hence, 
NVT is expected to occur in long AF rings with large 
anisotropy. While QTM can be probed by macroscopic 
magnetization measurements, NVT is elusive and a 
direct experimental demonstration of its occurrence is 
lacking. Neel vector tunneling has been recently shown 
to occur in a Feig ring Q, but only indirectly through 
powder INS and bulk thermodynamic measurements. 
The possible occurrence of the NVT scenario in the 
low-T dynamics of AF rings can be assessed by extract- 
ing dynamical correlations functions with the method 
presented in the paper. Indeed, in the case of NVT 
the two lowest eigenstates are close to the symmetric 
and antisymmetric superpositions of the two Neel states 
tit •■■) an d | tti •■■)■ This implies that all the spins 
are rigidly locked into a two-sublattice pattern. If a 
local perturbing field is applied along z on a single spin 
as discussed before, the expectation values of all the 
spins jointly increase in modulus (with an alternating 
pattern) with a practically instantaneous propagation 
of the perturbation. This behavior, shown in the lower 
panels of Fig. |U is completely different from the one we 
have experimentally found in Cr s (upper panels of Fig. 
2|). Hence, these results demonstrate that the low-T 
dynamics of Cr§ is not characterized by NVT. Figure 
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[4] shows that even for \D/J\ < 1 (but larger than in 
the case of Cr 8 ) NVT dynamics occurs. The present 
method could be applied to Fei§ to deeply investigate its 
spin dynamics and directly show that it displays NVT, 
circumventing also the fact that the huge dimension of 
its Hilbert space prevents exact calculations. 



V. CONCLUSIONS AND PERSPECTIVES 

By exploiting the Cr§ prototype AF ring as a bench- 
mark, we have shown that the quantum spin dynamics 
of MNMs can be fully unraveled by four-dimensional 
inelastic neutron scattering. The possibility to extract 
all two-spin dynamical correlation functions opens 
remarkable perspectives in the understanding of crucial 
but elusive quantum phenomena in several classes of 
magnetic molecules. In particular, supramolecular com- 
plexes containing linked molecular nanomagnets have 
been demonstrated to display entanglement [H, EH and 
to be excellent candidates for implementing quantum- 
computation Q and quantum-simulation algorithms 
@. The present method could be applied to reach a 
profound knowledge of the structure of the eigenstates 
of these complexes and of the entanglement between 
different molecular qubits. 

Another emerging field is that of frustrated MNMs, 
whose magnetic correlations and spin dynamics are 
particularly complex. In fact, their energy spectrum has 
a structure with large degeneracies which cannot be ra- 
tionalized in terms of simple paradigms (e.g., giant spins 
or rotational bands). However, a deep understanding 
of the spin dynamics of these systems is often hindered 
by the very large dimension of the associated Hilbert 
space. The approach described here can be exploited to 
determine the full pattern of dynamical spin correlations 
in some of these systems. 
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The original procedure for preparation of compound 
Cr 8 F 8 [(0 2 CC(CH3)3]i6[(CH 3 )2CO] 2 (1) reported in 
Gerbeleu, N. V. et al. Dokl. Akad. Nauk. SSSR 313, 
1459, (1990), and then the modified version published in 
[24l | produce the perdeuterated version of 1 in a very low 
yield (probably caused by the isotope effects on acidity of 
deutcratcd pivalic acid) and a new synthetic procedure 
for the deuterated compound CrgFg [(02CC(CD3)3]i6 
(2)was developed and is reported below. 
Heating in a Teflon flask at 150 °C for 14h while stirring 
of deuterated pivalic acid, ( (CD 3 ) 3 CC02H , 11.0 g, 
98.94 mmol), hexamethylenetetramine ( C6H11N4, 0.6 
g, 4.28 mmol) and chromium(III) fluoride tetrahydrate 
(CrF3 .4H2O, 3.0 g, 16.57 mmol) resulted in a green 
microcrystalline product. After this, the flask was cooled 
to room temperature, and acetonitrile (20ml ) was added 
to complete the precipitation. The product was collected 
by filtration, washed with acetone, dried, and then was 
extracted with hexane (50 mL), and finally purified on a 
silica gel column using toluene as the eluent. The first 
fraction from the column was collected and the solution 
was evaporated to dryness. Yield: 3.15 g (65 %). 
Elemental analysis calculated (%) for CsoCrgF8Di44032 
: Cr 17.85, C 41.22 ; found: Cr 17.80, C 40.99. ES-MS 
(sample dissolved in THF; spectra run in CH3OH) : 
2352+ [M+Na]+ , 2384+ [M+CH30H+Na]+. 
Large crystals of Crs were obtained by very slow evapo- 
ration of hexane solution of 2 at ambient condition. The 
crystal structure is tetragonal (space group P4212) with 
a '= b = 20.093(2)1 and c = 16.801(2)1. 

INS experiments were performed on the IN5 time-of- 
flight inelastic neutron spectrometer [l9| at the high- flux 
reactor of the Institute Laue-Langevin. The IN5 instru- 
ment has a 30m 2 position sensitive detector divided in 
10 5 pixels, covering 147 degrees of azimuthal angle and 
±20 degrees out-of plane. A 0.24 g Cr 8 single crystal 
was aligned as in Fig. SI (Supplementary Information) 
and measurements were taken by rotating the crystal 
(in steps of 1 degree) about the vertical z' axis. Incident 
neutron wavelength of 5.0 and 3.ll were selected to 
probe measured excitations with energy resolutions, at 
zero energy transfers, of 0.079 and 0.38 mcV respectively. 
The uniformity of the detector sensitivity was ensured 
by measurement of a standard vanadium sample. 
Measurements for different rotations an gles were com- 
bined using the HORACE analysis suitej2~L| . S(Q',uj) 
was rotated by 37.8 degrees about the Q y > axis to analyze 
the data in the (x, y, z) molecular reference frame. The 
two molecules in the unit cell (Fig. SI) were included in 
the calculation of the scattering function and fitting of 
Eq.([2]) to the Q dependence of the observed transitions 
was performed using HORACE. Some compensation 
for inhomogeneity in the scattering intensity due to 
self-shielding is performed, based on the variation of 
incoherent scatting intensity with sample orientation. 
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VI. SUPPLEMENTARY FIGURES 




Y 



FIG. 1: Orientation of the Crg single crystal for the IN5 ex- 
periment. (x',y') defines the laboratory horizontal scattering 
plane and (x, y) is the plane parallel to the ring. The dotted 
line is the unit cell. The crystal was mounted with the [1 
1] direction approximately perpendicular to the scattering 
plane. The angle between z, perpendicular to the ring plane, 
and the laboratory z' direction is 37.8 degrees. 




FIG. 2: Illustration of the equal-time correlations {sz(d)s z (l)) 
deduced from Table I of the main text and mapped onto the 
length of the cones. The red cone corresponds to d — 1. 




FIG. 3: Constant energy plots of the neutron scattering in- 
tensity for the two resolved excitations in the p=l peak (see 
the inset in Fig. 1 of the main text). As in Fig. 3 of the main 
text, we show the dependency of intensity on two wavevector 
components Q x — Q y lying in the ring's x — y plane, integrated 
over the full Q z data range, a (b): Data from the higher-E 
(lower-E) subpeak measured with a 5 Aincident neutron wave- 
length. The positions of the maxima are the same for both 
peaks but the relative intensities of the various maxima are 
different. In particular, the two Q x = maxima are the most 
intense for the 0.7 meV transition but are the weakest for the 
0.9 meV one. c-d: fits to Eq. 2 (coefficients are reported in 
Table II in the main text). 



